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It is shown here that surfaces of high-temperature superconductors are covered by dipole layers. 
The charge density modulation is induced by the local suppression of the gap function at the surface. 
This effect is studied in the framework of the Ginzburg-Landau theory and crucially depends on the 
appropriate boundary conditions. Those are derived from Gor'kov's equations for a d-wave pairing 
symmetry. Within this framework the structure of the surface dipole layer is determined. The 
contribution of this charging to a lens-effect of superconducting films with holes, which has been 
studied in recent experiments, is discussed. 

PACS numbers: 73.30.-)-y, 74.50.+r, 61.16.Bg 
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I. INTRODUCTION 



If superconductors were absolutely "perfect" conduc- 
tors, they should screen charges and would be completely 
free of internal electric fields. While this is true on macro- 
scopic length scales, it certainly can not be expected 
on atomic length scales. Some recent attention focused 
on the surprising fact that such electrostatic effects may 
even occur on scales of the correlation length, which typ- 
ically is considerably larger than an atomic length. 

Electrons in the superconductor are equilibrated if 
their efeciroc/iemjca/ potential is spatially constant. How- 
ever, the formation of a superconducting state is in gen- 
eral accompanied by a change of the chemical potential 
of electrons. In a spatially inhomogeneous situation the 
modulation of the chemical potential induces an accumu- 
lation of electric charge density 0-0] such that the result- 
ing electrochemical potential is constant. These charging 
effects strongly depend on the ratio of the gap A to the 
Fermi energy ep and therefore are more strongly pro- 
nounced in high-temperature superconductors (HTSC) 
than in conventional superconductors. 

Several aspects of charging effects have been examined 
since the availability of the HTSC: an anomalous temper- 
ature dependence of the work function [Bj , charge redis- 
tribution effects within the layered structure of HTSC, 
and charging of vortices (2[g] . As a direct probe for the 
latter effect Blatter et al. suggested [|| to observe the 
electric stray field near the surface by atomic force mi- 
croscopy. Unfortunately, at present the expected effects 
are beyond the resolution of this experimental technique. 

In order to investigate alternative possibilities to ob- 
serve such charging effects, we examine surfaces of HTSC 
in a vacuum. Our motivation is twofold, (i) The higher 
dimensionality of a surface compared to a vortex line can 
be expected to lead to the accumulation of much a larger 
charge quantities. Even if this does not necessarily lead 
to much higher electric field amplitudes, the field will be 
extended over a much larger region, (ii) In a recent ex- 
periment PI the influence of a thin superconducting film 
on an electron beam penetrating a hole in the material 



has been examined. At T^ a change of beam intensity be- 
hind the hole was observed, i.e. the hole effectively acted 
as a lens. In principle, charging induced by the suppres- 
sion of the order parameter at the surface near the hole 
could provide an explanation of this electro-optical ef- 
fect. This mechanism will be studied quantitatively in 
this work. 

During the last few years increasing evidence has been 
found for d-wave pairing instead of a conventional s-wave 
pairing in HTSC [g[. For this reason we take specifically 
account of a d-wave symmetry. In the case of a vor- 
tex line the gap vanishes for topological reasons in the 
vortex center for d-wave pairing as well as for s-wave 
pairing [Q. The structure of the vortex core does not 
feel the underlying symmetry of the order parameter, at 
least in the vicinity of Tc. However, at a superconduc- 
tor/insulator interface the strength of the suppression of 
A and eventually also the charging does crucially depend 
on the symmetry. 

In Section O we give a derivation of the gap profile 
near the interface for general singlet pairing on the basis 
of the BCS theory. In particular, we formulate appro- 
priate boundary conditions for the order parameter in 
a phenomenological Ginzburg-Landau description. The 
resulting charging effects at an infinite plane surface are 
calculated in Section H. The electrostatic calculation of 
the lens stren gth of a superconducting film with a hole 
follows in Sec. IV . Section ^ concludes with a discussion. 



II. 



BOUNDARY CONDITIONS FOR d-WAVE 
ORDER PARAMETER 



A spatially varying order parameter profile can be 
described in the framework of the phenomenological 
Ginzburg-Landau (GL) theory. The physics near the sur- 
face crucially depends on the imposed boundary condi- 
tions. As shown by Gor'kov, this phenomenological the- 
ory can be derived from the microscopic Bardeen-Cooper- 
Schrieffer theory [gj . We essentially follow this approach 
in order to determine the boundary conditions appropri- 
ate for HTSC. 



The fact that the tetragonal HTSC materials 
YBa2Cu307_5 and La2-a;Sr^Cu04 belong to the class 
of unconventional d-wave superconductors appears has 
been reliably established [|| . The term "unconventional" 
means that the spatial symmetry of the superconducting 
order parameter Ak,Q/3 ^ (akaO-k^) is lower than in the 
normal state |Q]. 

More specifically evidence favors a (ia;2_j,2 -symmetry. 
In this case the order parameter can be written in the 
form 



(a) 
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Ak,a/3(r) = (icry)a/5?/'(k)A(r) 



(1) 



(b) 



with the Pauli matrix Uy. The momentum dependence 
at the Fermi surface is described by the following nor- 
malized basis function of the irreducible representation 
Big of the tetragonal group D^h'- 



V'(k) 



(fc| 



rCv 



(2) 



where X and Y are internal axes of the crystal, and 

k = k/|k|. 

In the vicinity of T^ the gap profile A(r) near the su- 
perconductor/insulator interface is determined by the so- 
lution of the GL equation 
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Explicit values of the correlation lengths C||(_l)(2^) along 
different axes and the gap satur ation amp litude Aq are 
given in Appendices [Eqs. (AlC) and (B4)]. 



To obtain the boundary condition for the GL equation 
(0) microscopically, we start by noting that the Gor'kov 
equations take the form of a linearized integral equation 
near T, il: 



A(ri)= I d\2 K{vuV2)A{v2). 



(4) 



The kernel K can be calculated in the quasiclassical ap- 
proximation, making use of the method of classical tra- 
jectories B,[lO||. Depending on the roughness of the sur- 
face, different kinds of the reflection of electrons, usually 
referred to as diffusive and specular, are to be considered 
(see Fig. |). 



(c) 




FIG. 1. Electron paths of the quasiclassical approximation: 
(a) straight paths contribute to the kernel independent of 
the surface, (b) for specular reflecting surfaces one additional 
trajectory contributes (c) for difi'usively reflecting surfaces a 
whole group of additional trajectories contributes. 

We consider the following geometry: The normal vec- 
tor n to the superconductor/insulator interface lies in the 
basis XY plane, the angle between n and the axis OX of 
the underlying tetragonal lattice being equal to $. The 
interface is assumed to be macroscopically flat, i.e. the 
roughness is restricted to scales smaller than the correla- 
tion length. We denote ^ = ^|| unless stated otherwise. 

Let us start with the case of diffusive reflection for 
microscopically rough surfaces. In the absence of an ex- 
ternal magnetic field the order parameter depends only 
on the distance x from the surface. It is convenient 
to introduce the dimensionless coordinate x = x/^q 
(fo = wf/SttTc is the coherence length), and the gap 
equation (Eh takes the form 



A(ii) = / 
Jo 



dx2 K{XI,X2)A{X2) 



(5) 



where 
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where V is the coupUng constant, and A^o is the electron 
density of states. The kernel K in (0) is composed of a 
bulk and a surface contribution. Both crucially depend 
on the angle $ through the functions Fq and F^: 



Fo (s) = i| (3 - Us^ + lOs^) cos^ 2$ 



-Fr(si,S2) = ^(l-3s?)(l-3s^)cos2 2$. 



(7) 



(8) 



Details of the derivation of these expressions can be found 
in Appendix O. 

In order to obtain the boundary condition we have to 
evaluate the linearized gap equation (m with the kernel 
(^. We will do this at the critical temperature Tc, where 
the BCS condition 



+00 



dx2 Kq(xy - X2) = 1, 



(9) 



holds. Here Kq is the bulk part of (||). In the region 
1 <^ X <^ ^{T)/^o the surface contribution of the ker- 
nel is small compared to the bulk contribution and the 



non-linearity of the gap equation still can be neglected. 
Therefore the gap equation is solved by a linear function, 
A{x) = A(0)(H-a;/6), so that the effective boundary con- 
dition for the order parameter in the GL region takes the 
form 



dA 

dx 



x=0 



u 

b 



(10) 



5=0 



The "extrapolation length" b acquires from the kernel a 
dependence on the orientation of the surface with respect 
to the underlying crystal lattice. 

We evaluate the parameter b using the variational ap- 
proach of Ref. ||l^. It is convenient to introduce the 
function q(x) by A(i) = C{x + q{x)), from which the 
parameter b follows according to 6 = lim q{x). The gap 

X — >oo 

equation can now be rewritten as 

9(2;i) = o^(^i) + / dx2 K{xi,X2)qix2), (11) 



where 



F(xi 



dx2 X2K{xi,X2) - 2xi 
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Apart from a prefactor, the solution of Eq. (|11|) can be 
obtained by minimizing the functional 



n<i] = 



/p°° dx q{x) [q{x) - /„°° dx'Kjx, x')q{x')\ 
[/o°° '^^ (i{i)E{x)\ 



with respect to q. The minimum value is given by 

1 
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2j^ di E{i)q{i)' 



(12) 



(13) 



The extrapolation length can be related to this mini- 
mum value as follows: Eq. (|ll|) can be rewritten as an 
equation for Q{x) = q{S:) — &, which vanishes for x — > 00. 
Upon multiplying both sides of the equation for Q by 



xi, integrating over xi and taking into account (|13D, we 
finally obtain jl2| (see also ]lq|) an exact expression for 
the extrapolation length b — b^Q: 



Co 7C{3) J^ ds s^Fois) 
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+ dsi ds2 slslFr{si,S2) 
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ds s^Fo{s) 



1 



VNoJ^n 



(14) 



where C{x) is the Riemann zeta- function. 

Now we are able to apply the variational principle. 
Substituting in (n2h unity as a trial function, we have: 
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0-' mill 



27r^ /o ^'5 s^o(s) - /o c^si /o '^^2 siS2-FV(si, S2) 



49(2(3) 



Jg ds s2Fo(s) + Jg dsi /(J ds2 sis^i^r(si, S2) 



(15) 



As one can see from Appendix M, the correspon ding 
expressions for the specular case follow from Eqs. (|lj) 
and (^ by replacing F^{si,S2) with ^(5(si - S2)Fi.(s), 
where i^r(s) is given by 
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19s'') cos^ 2$ 
-s'*)sin2 2$. 



(16) 



Then the results of Ref. ||ll[ are recovered. 

The extrapolation lengths for the diffusive and spec- 
ular cases, which result from the substitution of expres- 
sions (0), (p^, and (||) in (|l5|), are plotted in Fig. g as a 
function of the angle $. 



rfA 
dx 



= 0. 



(18) 



x=0 



In the consideration given above, the normal vector n 
was assumed to lie in the basis plane of a tetragonal crys- 
tal. If n is directed along the fourth order axis, then the 
functions Fo and F^ take the following form: 



Fo(s)=F,(s) 
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(l-s^)^ 



Fr(si,S2)=0 



(19) 
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FIG. 2. Extrapolation length b as a function of the angle 
$ for the diffusive (solid line) and the specular (dashed line) 



The profile of the order parameter near the surface 
is determined by the ratio of the extrapolation length b 
and the correlation length f (T). In the diffusive case the 
value of the extrapolation length b just slightly oscillates 
as a function of $. For all orientations of the surface it is 
of the order of the coherence length ^0 • In the GL region 
it is therefore much smaller than the characteristic scale 
of the order parameter, the correlation length ^(T), and 
the boundary condition for the GL equation effectively 
becomes 



AL 



0. 



(17) 



In the specular case the ratio 6($)/Co strongly oscil- 
lates from to 00. However, in the GL region 6(<I>) ^ 
£,{T) for most orientations, and effectively A 1 2,^0 — 
again. Only in a narrow range of orientations 5 3> ^, 
where the boundary condition is effectively 



[see (H), ( |b1i1) and (^13|)]. As follows from (|l|), the 
exact boundary condition along z-direction for specular 
reflection is (cf. Refs. |Jl|]) 



rfA 

dz 



= 0. 



(21) 



z=0 



In the diffusive case the corresponding gap equation can, 
in principle, be solved exactly since it has a Wiener-Hopf 
form due to (pO|). However, we shall not proceed in this 
rather cumbersome way but instead use the variational 
method again, which gives b k, 0.46^o. 

Thus, the final conclusion is the following: We are able 
to use the boundary condition (|l^) for all orientations of 
a diffusively reflecting surface with a normal vector lying 
in the XY plane. For a specularly reflecting surface this 
condition holds for almost all orientations, except from 
very narrow angular regions near $ = 0, ^ etc, where we 
should use (^8|) instead. The width of these regions is of 
the order of ^o/^{T), which is negligibly small near T^. 

The profile of the order parameter at the surface is de- 
termined by the solution of the GL equation (^ subject 
to the boundary condition (ITol). One obtains 



A(x) = Ao tanh[(a; - xo)/\/2^(T)] 



with 



2^0 



-\/2^(T)arctanh 



e(T) 

%/26 



e{T) 

2&2 



(22) 



(23) 



In the case of a diffusively reflecting boundary we have 
obtained 6 ~ ^o- Thus, in this case it is reasonable to use 
the estimate a;o ~ near T^. 



III. SURFACE CHARGING 

To obtain the charge modulation near an infinite plane 
surface we follow the calculations for the charging of a 
flux line by Blatter et al. p] . The spatial variation in the 
order parameter induces a modulation of th e ch emical 
potential ^. As derived in Appendix ^ [Eq. ( Al3| )], this 
generates a variation of the local charge density 



Qc^t{x) 



^eN',\n^[l^l 



-4e7VoAg In 



Tc 



^\x) 



exp 



%/2a 



(24) 



where N^ is derivative of the density of states at the 
Fermi level, Ec is the energy cutoff of the electron in- 
fraction, Aq denotes the bulk value of the gap function, 
and A^o denotes the density of states at the Fermi level. 
We consider Qcyit{x) as "external" charge density since 
screening has not yet been taken into account. 

Metallic screening can be included within a Thomas- 
Fermi approximation by introducing the spatially con- 
stant electrochemical potential fj, + ecj) instead of /i. The 
electrostatic potential is then determined in linear or- 
der by the one-dimensional screened Poisson equation: 






1 



(x) = -i-KQc^tix). 



(25) 



The Thomas-Fermi screening length is given by Atf = 
(87re2Ao)-i/2. 

Overall charge neutrality in the superconducting half- 
space a; > requires a solution that meets the boundary 
condition ^'(a; = 0) = 0. In the limit Atf ^ £, this 
requirement is fulfilled by the solution 



(a;)«-167reA^A2^1n 



^/2\ 



tfg 



-x/\i 






-V2x/^ 



(26) 



(28) 



per unit surface area. Here we have used that h\{ec/Tc) 
is usually of order of unity for high-Tc superconductors. 
The corresponding dipole moment g per unit area of the 
surface is given by 



9 



xg{x)dx : 



-^* 



(29) 




FIG. 3. Modulation of the gap A{x) and the charge density 
p{x) at an infinite plane surface as a function of the normal 
distance x from the surface (arbitrary units). 



IV. LENS EFFECT 




layer of dipole s 



The corresponding screened charge density is g{x) = 
—(J)"{x)/Ati. In t he ca se of d-wave pairing one gets, with 
the specific gap ( |A10D and correlation length (B4) for a 
three-dimensional parabolic band, the explicit expression 



g{x) 



1 eoB 
5^ 



In 



V2 



A- 



c/Atf „ ^g-V23;/« 



TF 



i 



(27) 



with the Bohr atom radius ob. 

Due to screening this charge distribution forms a dou- 
ble layer of opposite charge, see Fig. ^. The outer layer 
with a thickness of the order of Atf ln(^/V2ATF) con- 
tributes a total charge 



FIG. 4. Superconducting film with circular hole inside act- 
ing as lens for a penetrating electron beam. 



So far we have derived that the surface of a high-Tc 
superconductor can induce a charge modulation. Now we 
comment on a possible experimental observation of this 
effect. One setup, an "electrooptical lens", is shown in 
Fig. ^ and can be considered as idealization of a recent 
experiment Q. Motivated by the observed lens effect, 
which sets in when the material becomes superconduct- 
ing, we will examine to what extent it can be ascribed to 
a surface charging according to the presented mechanism. 

Let us assume that the sample is a thin d-wave super- 
conducting monocrystalline film with a micro-size circu- 
lar hole inside. The surfaces of the film are assumed to be 



smooth and parallel to the XY crystal plane. Assuming 
specular reflection, the order parameter is not suppressed 
at these surfaces according to Eq. (|2l]). But the inner 
surface of the hole is normal to the XY plane leading 
to a strongly suppressed gap function at the boundary 
of the hole. Due to the preparation of the hole this sur- 
face should be rough leading to a decrease of the order 
parameter around the hole according to Eq. (17). Thus 
one obtains a charge accumulation around the hole only 
|l4| . Since the real charge distribution consists of two 
layers of positive and negative charge one can think of 
this distribution as a ring of dipoles oriented normal to 
the boundary of the hole, see Fig. |^. The corresponding 
electrostatic field should be observable in different ways. 
One possible way to detect the ring of dipoles is to send a 
very slow and weak beam of electrons through the hole. 
The intensity in the center of the beam image behind 
the hole then changes when the sample is cooled below 
Tc due to the (de)focusing caused by the accumulated 
dipoles. 

In the actual experiment by Kriebel et al. |5[| a simi- 
lar geometry was examined by transmission electron mi- 
croscopy. They have observed a very small change in 
the image intensity of the electron beam at the super- 
conducting transition. However, the actual sample was 
polycrystalline and can not be expected to lead to such 
a well-defined field distribution as in our ideal situation 
described above. Nevertheless, to obtain the order of 
magnitude of this lens effect, it is sufficient to consider 
the simplified setup. 

The calculation of the lens effect amounts to solving 
the electrostatic potential problem for the ring of dipoles. 
From the electric field we calculate the optical effect in 
an eikonal approximation for the electrons. 

Since we can assume that the radius of the hole is much 
larger than the correlation length ^, the curvature of the 
boundary can be neglected. Therefore the accumulated 
charge can be approximated by the charge density ( |27| ) 
obtained for a plane surface. With the (constant) dipole 
moment g per unit area from Eq. (E3) the electric po- 
tential outside the superconducting film is given by 



«^£ 



n- V 



1 



dS' 



(30) 



with the surface normal vector n lying in the XY plane 
and the surface element dS' of the ring. For a cylinder 
with radius R and height d -^ R this leads to 



(j){p, z) 



gd 

An 



dip 



Rp cos (f ~ R^ 



[p2 + z2 + i?2 _ 2Rp COS ip] 



3/2 • 



(31) 



This integral can be expressed in terms of complete el- 
liptic integrals. But we are only interested in the radial 
electric field component averaged along the z direction to 
obtain the total intensity change. Let us assume that an 
electron moves initially at a radial distance pi from the 
center in the negative z-direction with vanishing radial 



velocity. To minimize the influence of the source of the 
electron on the effect one should assume that the source 
is asymptotically far away. 

For this radial distance pi from the z-axis the averaged 
value of the deflecting field for — L < z < 0, defined by 



Ep{pi) = - / dzEp{p,,z), 



(32) 



can be calculated as follows: Consider a virtual cylinder 
of radius pi around the z-axis with its ends at z = and 
z ^ —L, respectively. For symmetry reasons the flux of 
the electric field across the abutting face at z = van- 
ishes. If the lower face of the cylinder is far away from the 
dipoles, i.e. L ^ R, the electric potential at the lower 
face can be approximated by an asymptotical expansion 
of ( pl| ) in z. Then the angular integration is trivial and 
one obtains for the z-component of the asymptotic elec- 
tric field 



E.{p,z) 



3gdR^ 2z2 - lOp^ - 5R^ 



+ 0(z-8). (33) 



Because there is no charge located inside the virtual 
cylinder, the total fiux of the electric field across the sur- 
face of the cylinder vanishes. Using the asymptotic value 
( P3| ) to calculate the flux across the face at z = — L one 
obtains the averaged radial field for z < 



EpiPt 



-gdpiR 



, 5{pI+R^)-2L^ 



(34) 



For symmetry reasons the same expression holds for an 
average of the electric field in the upper half space. We 
take L ^ oo as starting position of the electrons. There- 
fore the electric field in the upper half space does not 
contribute to the total deflection of the electron, as can 
be seen easily from (34). Since the whole change in the 



intensity comes then from the path of the electron behind 
the film, the optimal position z = —L of the detector has 
to be obtained by a maximization of the averaged ra- 
dial field (p3) as a function of L. For electrons near by 
the center of the hole {pi « 0) the optimal distance is 
L = (\/T4/2)i?. For this rather small value of L it can be 
checked numerically that the asymptotic approximation 
used for the average procedure underestimates Ep only 
by a factor « 2 if pi <^ R. Using this optimized value, 
the electron will be detected at z = — L with a radial 
distance Pf{pi) from the center given by 



PfiPi) = P»(l-ao + a2pD 



with 



ao ~ 0.016 



a2 « 0.04- 



gd 



m Rvq ' 



gd 
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(35) 

(36) 

(37) 



where vq is averaged velocity of the electron with respect 
to the z direction. The relative enhancement of the in- 
tensity at the center is then given by 



dpf 



1 



p/=o 



1 - ao 



Using the density of the dipole moment (|2 
for the deflection parameter 



(38) 



we obtain 
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e^flB 
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IQ-^eV d 
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(39) 



where we have used ^ ~ lOA as a typical value for high- 
Tc superconductors in the last expression. The spatial 
variation of the relative intensity I{pf) — dpi/dpf deter- 
mined by (35) is shown in Fig. ra. Strictly, the shown 



numerical result is valid only for distances pf <$^ R due 
to the truncation of higher order terms of pi in Eq. (35). 




0.4 0.6 

Pf/R 
FIG. 5. Variation of the relative intensity I{pf) of the elec- 
tron beam image across the hole behind the superconducting 
film. 



V. DISCUSSION 

In Sec. |l^ we have presented an explicit microscopic 
calculation of the extrapolation length in d-wave super- 
conductors. Within the quasiclassical approximation we 
found 6 ~ <^o, which is of the order of 12 A in YBCO. 
Within the same approximation microscopic calculations 
||15| for s-wave superconductors give b = cxd, i.e. no sup- 
pression of the order parameter at the surface. In the 
s-wave case a finite value of 



b ^ eo/a 



(40) 



can be found only if one goes beyond the quasiclassical 
approximation, taking special care of the finite thickness 
of the metal- vacuum transition layer ||l^,|l5[ . From this 
expression one can expect small 6, i.e. a notable suppres- 
sion of the condensate near the surface, if i^o is small (i.e. 
of the order of a). This is the case for HTSC as has been 
pointed out by Deutscher and Miiller [nTl. However, Eq. 



( pt]| ) is based on approximations which are valid only as 
long as 5 3> ^0- For d-wave superconductivity the situa- 
tion is very different: as our calculation shows, the most 
important reason for the reduction of b is the d-wave sym- 
metry itself. In this case the effect is so pronounced that 
it can be obtained within the quasiclassical approxima- 
tion. 

It should be emphasized that in this article we assume 
the superconducting order parameter to have the pure 
d3.2_„2-symmetry [see Eq. (0)], which immediately fol- 



lows from the chosen form of the pairing potential (A3) 
containing the spherical harmonics corresponding only 
to the irreducible representation Big of the tetragonal 
group. In general, one could allow for the possibility of a 
small admixture of the s-wave component. As shown in 
Ref . p8| , such an admixture appears whenever one deals 
with a spatially non-uniform distribution of the order 
parameter, e.g. near surfaces, twinning planes, columnar 
defects etc. and is brought about by the presence of the 
mixed-gradient terms in the GL expansion of the free en- 
ergy. In our work wc neglect this effect due to the results 
of Ref. |l^ : (i) The induced s-wavc component is always 
much smaller (by one order of magnitude at least) than 
the d-wave one. This holds in particular at temperatures 
close to Tc, since the s-wave component vanishes faster 
with (Tc — T) than the d-wave component, (ii) The s- 
wave component vanishes at the surface together with 
the d-wave component, so that the modulation of p(x) 
shown in Fig. ^ remains unaffected. 

The fact, that the d-wave order parameter is 
strongly suppressed at diffusively reflecting superconduc- 
tor/insulator interfaces, could explain an anomalously 
weak temperature dependence of the Josephson critical 
current jc ~ (Tc — T)^ near Tc, as observed in Ref. 
|19[ . Indeed, the critical current is proportional to the 
square of the order parameter at the surface: jc ~ A^ g. 
However, as it follows from our results [see Eq. (p2)]: 
Ajs -- (Ve(T))Ao ~ {Co/aT))Aa. Hence 



AliT) 



?o 



e{T) 



(Tc - Tf 



(41) 



for all orientations, instead of jc ~ Tc — T, as expected for 
a conventional SIS junction. The same result has been 
obtained in pO| for distinguished orientations of a spec- 
ularly reflecting interface (see also the discussion in the 
end of Sec. ||). 

In the present analysis we have neglected possible 
anisotropy of the Fermi surface, which is characteristic 
for HTSC. In terms of the ratio of effective masses in a 
Ginzburg-Landau theory this anisotropy can be of the 
order mj_/m\\ ~ 100. The generalization of our analy- 
sis to anisotropic Fermi surfaces is straightforward but 
rapidly leads to cumbersome analytic expressions. Such 



a generalization has been performed by Shapoval |21| for 
the case of an anisotropic s-wave pairing. For d-wave 
pairing we expect that this generalization leads only to 



quantitative changes which only weakly affect the bound- 
ary condition b <C ^{T). Although N[) and Atf ~ N^^^'^ 
drastically depend on the mass anisotropy. In contrast 
the combination -/VpA^p is insensitive. Therefore only the 
spatial distribution of charge g will be affected, but not 
the total charge Q or the dipole moment g. 

The amplitude of the electrooptical lens effect, Eq. 
(p9|), crucially depends on the beam energy of the elec- 
trons. In the experiment of Ref. [^ a beam energy of 
E\cin — 120keV was used, which leads to ao ~ 10~^ at the 
optimum distance L ^ R behind the lens. However, the 
actual picture was generated much further away, where 
the intensity rapidly decreases according to Eq. (p4|). In 
addition, the sample had a much more irregular shape 
than in our idealized geometry. We expect this irregular- 
ity to reduce the electrooptical effect of the dipole layer. 
Therefore we believe that surface charging has to be ex- 
cluded as the origin of the lens effect with ao ^ 10^^ 
observed in this particular experiment. 

At present, the actual origin for the observed effect 
can not be deduced unambiguously from the experiment. 
One contribution could be due to a charge accumulation 
at the surface of the superconductor which is directly 
hit by the electron beam. The amount of accumulated 
charge then strongly depends on the conductivity of the 
sample which drastically changes at Tc- Since in the ex- 
periment polycrystalline samples have been used, there 
are small magnetic fields induced by the Josephson cur- 
rents between the grains. But these fields should give 
an even weaker deflection of the electrons since the local 
field sources do not sum up constructively. 

In order to clarify this situation it would be desirable 
to use lower electron energies where the the electroopti- 
cal effect should become stronger. For i^kin = lOOeV as 
used in LEED, ao ^ 10~^ can be achieved due to surface 
charging at the optimum distance. In addition, the use of 
a single crystalline sample with a more regular shape of 
the hole would be desirable as well as a spatial resolution 
of the density profile. 

Finally, to detect the charging effect, other experimen- 
tal techniques might be more promising than the lens 
effect examined above. For example, the local stray field 
outside the superconductor may be probed by atomic 
force microscopy, as already suggested by Blatter et al. 
M for vortices. 
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APPENDIX A: d-WAVE PAIRING 

We briefly point out the most important changes for 
pairing in contrast to conventional s-wave pairing 
in a weak-coupling BCS model. A discussion of arbitrary 
symmetry has been summarized in Ref. [BJ. For calcu- 
lational convenience we ignore mass anisotropy and use 
a dispersion relation ek — k^ /2m — /i of single electron 
states. They have a density of states (per spin projection) 



iV. := 
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27r2 



Irrespective of specific pairing mechanism, the gap 
function for singlet Cooper pairs with relative momen- 
tum k can be written as 

f (Pk' 
Ak,Q/3 = (icr;^)a/3V'(k)Ao = / V (k, k') (ak'aa_k';3)- 

(A2) 



by 



For a singlet state the potential can be approximated 



nkk') = /"^^^^)^'*^^')' l^k|,|ek'|<ec 
^ '• ' 1 0, otherwise 



(A3) 



with an energy cutoff ec of the interaction. For 
anisotropic masses the definition (0) of i}} was normal- 
ized such that (■02(k))o = 1 with an average 
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defined on the Fermi surface. 

In the BCS mean-field (MF) approximation the Hamil- 
tonian is diagonalized by the Bogoliubov transformation 
with parameters 



|wkP = (l + ek/^k)/2 
|«k|' = (l-ek/^k)/2 



Ei.= 



- . L2 



|Ao|2^2(k) 



(A5) 
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The self-consistency equation for the MF approximation 
reads 



l^V 



(fk 
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(A8) 



with the Fermi-Dirac distribution function f{E) = 1/[1 + 
exp(£'/T)]. The primed integral runs over k with |ek| < 
Ec only 

For Ao ^ Eq. ( |A8|) implicitly determines 

rc = 1.14eeexp(-l/T/7Vo). (A9) 

Due to the choice of the normalization of V', this relation 
does not depend explicitly on pairing symmetry or mass 
anisotropy. 



In the critical region the amphtude of the gap follows 
from Eq. (AS) after expansion in second order of Aq and 
linearizing in 1 — T jT^. 



|Ao(T)p = (V^V^(i~r/^c)r2. 



(AlO) 



In comparison to the isotropic s-wave case this result is 
modified by a factor (V''')o = 15/7. 

We now calculate for fixed electrochemical potential 
the change in electron density due to formation of a gap. 
In the superconducting state a single electron state k has 
occupation number 



nL = kkl' + (kkl' - Kn/(i?k) 

compared to the normal state with nj^ = 
change of the total electron density 



bn = 2 



(2^ 



(K - <) 
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/(ek). The 
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can be evaluated using the Sommerfeld approximation 
iVf « A^o + e^o- C)ne finds to lowest order in Ao, i.e. 
close to Tc, 



fc«iV^|Aopln-f, 



Tr 



(A13) 



where again all dependence on pairing symmetry and 
mass anisotropy is contained in Tc. Since the conden- 
sate density |Aop has a discontinuous slope near Tc, Eq. 



( Al3| ) implies that for fixed electrochemical potential the 
electron density has a discontinuous slope or, vice versa, 
the chemical potential for fixed electron density. This ef- 
fect and its consequences for the work function has been 
pointed out by van der Marel [0 . 



APPENDIX B: METHOD OF CLASSICAL 
TRAJECTORIES 

As shown in Ref. IqI, the kernel K in Eq. (0) is de- 
termined by the following expression, which is valid for 
arbitrary symmetry of the order parameter: 

/•oo 

i^(ri,r2) = TT^iVoT V / dt exp(-2|cc;„|t) 

„ "'0 



x(^*(ni)^(n2)),^, 



e—0, classical ' 



(Bl) 



where i > is the time of motion along the trajectory 
and LOn = (2n -I- V)t^T is a Matsubara frequency. The 
angular brackets denote the averaging over all possible 
classical trajectories r(i) of a particle, moving with the 
Fermi velocity vp^ which connect the points ri and r2 and 
satisfy the condition r(0) = ri (one has to take into ac- 
count both the trajectories, for which r(t) — r2, and the 
time-reversed trajectories, which give just a Hermitian- 
conjugated contribution to the kernel). The unit vectors 



ni and n2 denote the directions of velocity in the corre- 
sponding points. We shall keep below the general nota- 
tions for the basis functions (which are chosen to be real), 
so that our results are applicable for any one-dimensional 
spin-singlet order parameter. 

Let us first consider the contribution from straight tra- 
jectories. Th e expression in the angular brackets in the 
r.h.s. of (Bl) has the form: 
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J(r2-r(t))i/^2(n)-Hh.c., (B2) 



where r(t) = ri + vpvd is the only straight path from ri 
to r2 (see Fig. |^a, rii — v\.2 = n). Here, integration is 
over all directions of n. Hence 



ifo(ri-r2) = /fo(r) 



47reo r2 sinh(77^o) 



(B3) 



with ^0 = vf/'^t^Tc, and n = r/|r|. In the critical region 
the correlation lengths ^i along different directions are 
given by 



C'm=7, 



7C(3) 
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Due to the d-wave pairing, there is an anisotropy (even 
for the spherical Fermi surface) described by the factors 
Ix ^ lY ~ 9/7 and 7^ = 3/7. 

The equation (|^) follows from (|j) as a result of inte- 
gration over the differences 2/2 — 2/i and Z2 — zi. 



/•oo 

A(a;i) = / dx2 K{xi,X2)A{x2), 
Jo 



(B5) 



the kernel being given by the sum of two terms: 

K{xi,X2) = Ko{xi -X2) + K^[xi,X2), (B6) 

where the contributions Kq and K^ come respectively 
from the straight trajectories and those ones, which are 
reflected against the surface. 



Substituting (B2) in Eq. (Bl), we obtain: 



Ka{x) 
where 
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Here s — cos 6, 9 and if are the polar and 
aziniuthal angles respectively, so that n = 
(s, Vl — s^ cos ip, Vl — s^ sin ip). The polar axis is cho- 
sen along the normal vector to the surface (Fig. |^). 

Note that J^ dsFo{s) = 1 because of the normalization 
condition of the basis functions. 

The form of Kr{xi, X2) is determined by the reflection 
of electrons at the boundary. For the case of specularly 



reflecting boundary we have from (Bl) the following ex- 
pression: 

— 6(r2 - r(i))-0(ni)-0(n2) + h.c. 

The only possible reflected path is (see Fig. ||b): 



r(t) 



ri + vpnit 
VFn2{t - to) 



t <to = 
t>tQ , 



VpSi 



(B9) 



where to corresponds to the moment when 
the particle hits the sm'face and ni^2 = 
(±s, Vl — s^ cos (y9, Vl — s^ sin ip) (si < 0). Performing 
integrations we obtain: 



Kr{xi,X2) 



X exp 
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i(;{s,(p)ip{-s,ip). 
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The results (B8) and ( Bll ) obviously coincide with those 
obtained in Ref. fO] by another method. 

In the case of diffusive surface one needs to take into 
account all possible trajectories, reflected in all directions 
(see Fig. |^c): 
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ificctcd 



J ^ J dn2 P{ni,n2) S{r2 -m) 

X'tp{ni)^{n2) + h.c, 



where P(ni,n2) is the probability distribution for the 
particle to be reflected in the direction n2, provided it 
moved before along the direction rii. Assuming the re- 
flection is completely diffusive (i.e. the "memory" is com- 
pletely lost), we obtain that this probability density de- 
pends only on n2. As usual, it is convenient to choose it 
in the following manner: P(n2) = ^ cos 6*2 — ^S2- The 



corresponding trajectories are given by (B£), where ni 



and n2 are now independent unit vectors. After integra- 
tions, the resulting expression for K^ is |l2|: 



Kr{xi,X2) 



where 
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Using the specific angular depend ence (g) of the order 
parameter, the general expressions ( |B8| ) and ( |B13| ) lead 
to Eqs. (0) and (|). 

Note that if the contribution from the diffusively re- 
flected trajectories vanishes then the integral equation 
(^ can be solved exactly using the Wiener-Hopf method. 
This takes place for any spin-triplet order parameter, in 
particular for jp-w ave order parameter in the superfluid 
phases of ^He [|22| . 
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